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We investigate several models of coupled scalar fields that present discrete Z%, Z2 x Z2, Z3 and other 
symmetries. These models support topological domain wall solutions of the BPS and non-BPS 
type. The BPS solutions are stable, but the stability of the non-BPS solutions may depend on the 
parameters that specify the models. The BPS and non-BPS states give rise to bags, and also to 
three-junctions that may allow the presence of networks of topological defects. In particular, we 
show that the non-BPS defects of a specific model that engenders the Z3 symmetry give rise to a 
stable regular hexagonal network of domain walls. 

PACS numbers: 11.27.+d, 11.30.Er, 11.30. Pb 



I. INTRODUCTION 

Bogomol'nyi, and Prasad and Sommerfield jij, have 
presented a way to find topological solutions that sat- 
urate the lower bound in energy. These solutions are 
known as BPS states, and they are supposed to play some 
role in the development of modern particle physics, be- 
cause there are supersymmetric models where they are 
shorter multiplets, stable, and no continuum variation of 
parameters can modify these BPS states. In Refs. Q-^] 
one has investigated some models of coupled real scalar 
fields in bidimensional space-time. These models present 
peculiar features, such as the presence of the Bogomol'nyi 
bound, with the corresponding appearance of first-order 
differential equations that solve the equations of motion 
in the case of static configurations. These Refs. §-0] 
have found topological solitons in specific models by tak- 
ing advantage of the trial-orbit method first introduced 
by Rajaraman in Ref. ||. These investigations provide 
a concrete way of finding BPS states and suggest several 
other studies, in particular on the subject of defects in- 
side defects, as presented for instance in Refs. ||-|| and 
more recently in [PHl5f . The investigations of defects in- 
side defects presented in Ref. 0,1-2 14, lq] are natural 
extensions of the former papers [Q- |]7 

In supersymmetric models with chiral superfields the 
presence of discrete symmetry may produce BPS and 
non-BPS walls. The non-BPS walls are solutions of the 
equations of motion, and they do not necessarily saturate 
the lower bound in energy. But the BPS walls saturate 
the lower bound in energy, and lie in shorter multiplets, 
in a way such that they preserve the supersymmetry only 
partially |l6|-|T8|] . There are BPS states that preserve 1/2 
and 1 /4 of the supersymmetry fl6|-f23f , and in the present 
work we investigate several explicit examples. We start 
by first offering general considerations in Sec. 0. We ex- 
amine specific models in Sec. Ill, where we find several 



tigations concerning the BPS domain wall junctions in 
supersymmetric models. 

Our work continues in Sec. IV, where we investigate 
stability of several BPS and non-BPS solutions. This is 
done not only to show stability, but also to see how the 
bosonic matter behaves in the background of BPS and 
non-BPS defects. The behavior of the fermionic matter 
in the extended supersymmetric models can be investi- 
gated by following the lines of Refs. @,§7|,f||. In Secjv] 



VA 



we investigate two other issues. The first in Sec 
where we consider the entrapment of the other field 
in bag models p9|-|3l[|. This possibility may be imple- 
mented in the first, second and third models, and the 
bag may be of the BPS or non-BPS type, depending of 
the particular model u nder consideration. We examine 
the second issue in Sec. VB, where we investigate fusion 
of defects. Here we briefly review Ref. |2(J, which in- 
troduces investigations valid within the supersymmetric 
context, and after we offer another possibility of showing 
fusion of defects and stability of junctions, valid beyond 
the conte xt o f supersymmetry. 

In Sec. VI we show how one can generate planar net- 
works of BPS and non-BPS defects. There we show that 
the three-junction that appears in supersymmetric mod- 
els may be stable, but cannot generate a BPS network 
of defects. However, when we give up supersymmetry we 
find a model where domain wall junctions are stable, the 
three-junction strictly obey the triangle inequality and 
generate a stable network of defects. We summarize the 



main results in Sec. VII, where we end the present work. 



solutions, of the BPS and non-BPS type. See also the 
more recent work l23£9l , where one can find other inves- 



The subject of this work may find direct applications in 
diverse branches of physics. Besides the recent interest, 
put forward in Refs. pi|-p3[, we recall for instance that 
the Z3 group is the center of SU(3), and this may guide us 
toward applications involving strong interactions. Some 
examples can be found in Ref. j3^j3j|. Another possi- 
bility is related to the Z2 X Z3 symmetry, in an effort 
to entrap the ^-symmetric portion of the model inside 
the topological defect generated by the Z2 symmetry. 
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Other applications include issues related to cosmology 
|34| and magnetic materials |35|, and to pattern forma- 
tion in systems of condensed matter |36|. The investi- 
gations that we present are done in flat space-time, in 
the (d, 1) dimensional case. We use natural units, taking 
% = c = 1, and the metric tensor presents signature in 
the form (+, — , — ). In the present work we are mainly 
concerned with the cases d — 1,2, and 3. We specify the 
spatial dimension along with the interest of the subject 
under investigation. 



II. GENERAL CONSIDERATIONS 

Let us start with two real scalar fields <fi and \. We 
consider models defined via the Lagrangian density 



(1) 



This is valid in (d, 1) space-time dimensions. V = V((f>, x) 
is the potential. It has the general form 



v(<p, x ) 



1 



Wi 



-W 2 

2 x 



(2) 



W = W(<fr, x) is some smooth function of the fields 4> 
and x, and — (dW/d(j)) etc. This is the form of 
the (real bosonic portions of the) models we deal with in 
the present work. These models admit a supersymmetric 
extension, and in the supersymmetric theory the function 
W is the superpotential. 

Models defined via W(<fi, x) present several interesting 
properties . The equations of motion are 



d 2 4> 
W 



dV 







f?!* _ v 2 + — - 
dt 2 ^ dx 



(3) 
(4) 



For static field configurations living in the (1, 1) space- 
time, that is, for fields = 4>(x) and x = x( x ) they 
reduce to 



dx 2 
fx 
dx 2 



WfWw + w x w x4> 

WfW^ + W X W XX 



(5) 
(6) 



These equations are solved by first-order differential 
equations 



fx= W * 
dx * 



(7) 
(8) 



The energy of solutions that solve the first-order equa- 
tions are bounded to the value E l g = |AWiy|, where 



(9) 



Here (<f>i, X i) is the i-th vacuum state of the model under 
consideration. The reason for this is simple: the energy 
of static configurations is 



E 



dx 



dx 



dx 
dx 



wi + wi 



(10) 



It can be rewritten in the form 



dx 



dx 



+ 



dx_ 

dx 



We recall that 



dW 

dx 



d<t> 
dx 



dx 
dx 



(11) 



(12) 



to see that the energy gets to the bound Eb when the 
field configurations obey the first-order equations (^) and 
(H). This is the Bogomol'nyi bound, and the correspond- 
ing solutions are BPS solutions. The BPS solutions are 



linearly or classically stable - see Sec. IV 



The energy of BPS solutions can be also written as 
Em= I dx I I ^ ] + ( ^ ] I (13) 

(14) 

The above result shows that the gradient (g) and poten- 
tial (p) portions of the energy 




/oo 
dxg(x) = 
-oo 

p = / dx p(x) = 



dx 



dx — 

dx \ (W 2 



dx 
dx 



wi 



(15) 
(16) 



contribute evenly to the total energy of the solution. Al- 
though this is shared by all the BPS states, it is not 
peculiar to BPS states, since there are non-BPS states 
that also have energy evenly distributed in its potential 
and gradient portions - see below. 

In the case of two scalar fields, there are other possible 
arrangements to the energy (|l(]|). We can, for instance, 
consider the possibility 



E 



dx 



-j- — W x 

dx x 



dx 
dx 



(17) 



This expression is also minimized to obey Eqs. ( |13[ ) and 
(0) when one sets 



2 



dx x 
dx 



(18) 
(19) 



We can show that solutions to these equations (|8|) and 
( |l9| ) only solve the equations of motion when we impose 
on W((j), x) the extra condition 



XX 



(20) 



Unfortunately, however, this condition ( pp| ) factorizes 
W((f>, x) hito the sum W + {(f> + ) + W_(^>_) when one ro- 
tates the (4>,x) plane to (</>+, 0-), with 



(21) 



as firstly shown in Ref. [lq] . The condition (|20| ) turns 
the system into two decoupled systems. 

We can consider another possibility to rewrite the en- 
ergy @. We set 



E = 



1 



dx 



dd> 



dx 



-W4, 



(22) 



This expression is also minimized to obey (|lj) and Jl4| ) 
when we write 



dx x 
dx 



(23) 
(24) 



These equations (23) and (|24|) solve the equations of mo- 



tion when one imposes on W((f>, x) the extra condition 



w xx = o 



(25) 



This result is interesting. It shows that for W(4>, x) har- 
monic, we can solve the two (second-order differential) 
equations of motion (|^) and (||) with two sets of two 
first-order differential equations, the set of Bogomol'nyi 
equations (]?]) and (||), and another one, given by Eqs. ( p3| ) 
and (H). 

We make good use of this result by imposing that the 
second set of first-order equations (23) and (£4[) are still 
Bogomol'nyi equations, that appear via the presence of 
another function W . This new function is defined by 



(26) 
(27) 



Since W and W are smooth functions, they are neces- 
sarily harmonic. We then get the result that potentials 
defined by harmonic functions can always be written in 



terms of two functions, W and W . We use this to intro- 
duce the complex superpotential 



W = W + iW 



(28) 



The superpotential W = W(<p) now depends on the com- 
plex field tp = 4> + i x, and the harmonicity of both W 
and W appears via the Cauchy condition on W. This is 
the way we get from the investigations to the recent 
work in Refs. @-^3). 

We guide ourselves toward the topological solutions by 
introducing a topological current. This is the standard 
procedure, and we can introduce for instance the topo- 
logical current 



J% PS = e a Pd p W^,x) 



(29) 



It is defined in terms of W(</>, x), and the topological 
charge is directly identifyed with the energy of the topo- 
logical solution that satisfies the pair of first-order equa- 
tions (0) and (||). This definition identifies the topological 
charge with the energy of the BPS solution. It expresses 
the fact that the mass of a BPS state is completely deter- 
mined by its topological charge. This definition provides 
a single way to infer stability of the BPS states. This is a 
nice way of defining the topological current, but it is not 
complete, because it cannot see sectors of the non-BPS 
type. It is only effective for BPS sectors, and that is the 
reason for the label BPS assigned to it. However, since 
we need to identify every topological sector, we should 
use another topological current. We can consider, for 
instance 



ja = £ afi d 



(30) 



It obeys d a J a = 0, and it is also a vector in the (</>,%) 
plane. The topological current (|3(]) is the current that we 
shall use in this work. We use this definition to obtain, 
for static configurations 



JiJ a 



p 1 p 



2g(a0 



(31) 



where g(x) is the gradient energy density of the solution. 
This property can be used to infer stability of junctions, 
for configurations that present energy density evenly dis- 
tributed in their gradient and potential energy densities. 
This is so because in this case p l p = 2g(x) = g(x) + p(x) 
gives the full energy density, and may be used to show 
that the junction process occurs exothcrmically. This 
result was already introduced in Ref. [0, and offers a 
way of showing how the process of junctions of defects 
may occur exothermically in models defined outside the 
context of supersymmetry - see Sec. |V B| . 

Up to here the calculations were done in (1,1) space- 
time dimensions. However, we can easily immerse the 
system into (2, 1) or into (3, 1) space-time dimensions. 
In the first case, in d = 2 the solutions can be seen as 
domain ribbons, and the energy is now multiplied by the 
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length of the ribbons. In the second case, in d = 3 we get 
domain walls, and the energy is to be multiplied by the 
area of the wall. The planar case is necessary for intro- 
ducing junctions of domain ribbons [or of domain walls 
if one immerses the system in the (3,1) space-time] and 
their corresponding planar networks. There are several 
distinct possibilities, and in the sequel we specify some 
cases. We recall that we single out topological sectors 
using pairs of vacuum states. For BPS sectors the en- 
ergy is obtained via the superpotential, and we shall also 
use the notation = \W(<j>i,Xi) ~ W(c/)j,Xj)\i identi- 
fying the BPS sectors with tensions or energies of the 
topological defects. 



III. SOME SPECIFIC MODELS 

Our aim here is to introduce models of coupled scalar 
fields, to illustrate how one can find topological solu- 
tions of the BPS and non-BPS type. We investigate sev- 
eral models of coupled scalar fields, defined with smooth 
functions W(<j), x)- We organize the subject in the four 
subsections that follow, where we examine four different 
models. 



of minima of the potential. For X/fi > the system sup- 
ports several topological sectors. Their energies or ten- 
sions are given by, for solutions that connect the minima 
(-1,0) and (1,0) 



t 



ui-'S — g I'M 



(37) 



For solutions that connect the minima (±1,0) and 
(0,±y/X/n) we get 



t 



BPS 



|A| 



(38) 



For the non-BPS solution that connects the minima 
(0, y/X/fi) and (0, — a/A//x) we obtain 



t 



nBPS 




(39) 



There are BPS solutions in the first sector, which con- 
nects the two minima (±1,0). They are 



4>{x) = — tanh(A:r) 



(40) 
(41) 



A. The first model 



and 



We consider the superpotential 

W {1 \<i>, X ) =-X<i>+ l -Xct? +^X 2 (32) 

where A and /i are real parameters. This model was first 
investigated in Ref. Q], but here we bring new results. 
The equations of motion for static field configurations 
are given by 



= 2A 2 0(0 2 -1) + 2M 2 (2+^X 2 



fx 

dx 2 



2, ,3 



2A M (<r - i)x + Vrx + 



(33) 



(34) 



These equations are solved by configurations that also 
solve the pair of first-order differential equations 



dx 
dx 
dx 



X{4? -1) + mx 2 

2[i(j)x 



(35) 
(36) 



Solutions that solve these first-order equations are BPS 
solutions. 

This model has four absolute minima when X/fi > 0, 
and two minima when A//i < 0. For X/fi > the min- 
ima are (±1,0) and (0, ±y / X/fi). Topological solutions 
are solutions that connect adjacent minima. Thus, the 
asymptotic behavior of the fields is governed by the set 



<j)(x) = — tanh(2 fix) 

X(x) = ± 4 / — — 2 sech(2 u x) 
V A 4 



(42) 
(43) 



The first pair is of the one-component type, and is im- 
portant for investigating the presence of defects inside 
defects - see for instance Rcfs. OJi2|,|14|Jl5l. The second 



pair is of the two-component type, and may play a role 
in applications to condensed matter, as for instance in 
the case investigated in Refs. [||J37). 

The one-component and the two-component type of 
solutions also appear in magnetic materials J35| . There 
they are known as Ising and Bloch walls, respectively. 
The two-component Bloch wall has a peculiar behavior, 
that can be used to understand specific features of inter- 
faces between magnetic domains. The issue here is that 
if one looks at these solutions as vectors in the ((p, x) 
plane, for x spanning the interval (— oo, oo) they describe 
a straight line and an elliptical arc connection between 
the two minima (±1, 0), respectively, behaving very much 
like linearly and elliptically polarized light. For this rea- 
son, the Block walls can be seen as chiral interfaces be- 
tween domains with different magnetization, and this can 
be used to describe the wall motion in magnetic materi- 
als |38|. Furthermore, these walls can be charged and we 
may introduce charge by making the x field complex, or 
by coupling fermions to the system. 

For A//i < we have only two vacuum states. They are 
(±1, 0) and now we have a single topological sector, with 
the same energy and topological charge obtained before, 
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for the corresponding sector. We notice that the first pair 
of solutions (^) and ( fill ) solves the first-order equations 
even for X/fi < 0. This one-component solution is still of 
the BPS type for X/fi < 0, but now no topological defect 
can be nested inside it. The reason is that for A//i < 
the model presents no spontaneous symmetry breaking 
in the independent \ direction, and so it cannot support 
topological defect in this x direction anymore. 

For X/fi > 0, there are non-BPS solutions connecting 
(0, ^/X/Ji) and (0, — -\/A//i) by a straight line. To see this 
explicitly we consider (jj = in the specific system to get 
the equation of motion for the static field x = x( x ) 



xO) 



tanh(y Xfi x) 



(50) 



The pair <fi — and x m ( p>0| ) is the pair of non-BPS 
solutions that we have already found. Here, however, 
we found the non-BPS solution as an explicit solution 
of the new first-order differential equations in systems of 
coupled real scalar fields. Although this is a non-BPS 
solution, we have already shown in Sec. || that its energy 
is composed of equal portions of gradient (g) and po- 
tential (p) contributions. In terms of densities we have, 
explicitly 



dx 2 ^ 



The solutions are 



X(x) 




tanh(\/ X\xx 



(44) 



(45) 



gi(x) = pi(a;) = -A 2 sech 4 (^fxjlx^ 



(51) 



We can now understand this feature by realizing that 
although this pair of solutions is non-BPS, in the reduced 
system where (f> — > 0, the x ne ld solution is indeed a BPS 
solution. 



The corresponding non-BPS energy is (4/3)|A| \f~XjJi, as 
already presented in Eq. (|39|). 

This model decouples for A = /i, as we have already 
shown in Ref. |12|. In this case the symmetry Z2 x Z2 
changes to the Z± one, but here a rotation by it /A changes 
the model into two single-field models. We also note 
that for A = —2/j,, the model decouples into two single- 
field systems, showing that the limit /.i — ► —A/ 2 is also 
uninteresting, at least from the point of view of models 
of coupled scalar field. 

We expose other features of the model using We 
write 



W$ = 2A0, 



<t>x 



(46) 



We see that is harmonic when \i — > —A. This means 
that the model admits another pair of first-order equa- 
tions. This pair follows from (p3h and (p4|) 



B. The second model 

We consider another model, defined by the superpo- 
tential 



w^ 2 \4>,x) = X{4> 2 -i) x 



(52) 



This model was also considered in Ref. p8]. The equa- 
tions of motion are 



d 2 



d 2 



dt 2 dx 2 



2X 2 



1)0 + 4A 2 X 2 = O (53) 



and 



d 2 x d2 X.., 2 i2 n 



(54) 



For static fields the energy is minimized for configura- 
tions that obey the first-order equations 



dx 
dx 



= 

= A - Xcp 2 - ^x 2 



(47) 
(48) 



These equations also furnish solutions to the equations of 
motion, for fi — > —A. As we have already seen, for \x — > 
—A the model presents a single topological sector, already 
investigated. But we notice that W^ 1 ' is also harmonic 
in the limit (f> — > 0, independently of the parameters A 
and fi. We use this into the above Eqs. ( |47| ) and ( p8| ) to 
get (j> = and 



dx 



dx ^ \ /i 



A 



(49) 



We consider A^ > 0. In this case the above equation is 
easily solved to give 



-f = 2A</>X 
dx 

- - 1) 

dx 



(55) 
(56) 



The model presents two vacuum states, given by 
(±1,0). In contraposition with the former model, how- 
ever, this new model presents no BPS sector, because the 
function (|52|) obeys W^ 2 )(0 2 = l,x) = 0, and then gives 
|AW^ 2 )| = 0. This means that the first-order equations 
( p5| ) and (|5^) present no topological solutions connecting 
the two vacuum states. On the other hand, the vacuum 
states (±1,0) can be connected by topological but non- 
BPS solutions. The explicit form of the solutions is given 
by 



4>(x) = ±tanh(Ax) 
X (x) = 



(57) 
(58) 
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The energy or tension of the solution is 
tlsps = (4/3)| A| 



(59) 



This result reproduces the value presented by the similar 
one-component BPS solution (|l0|) and ( pll| ) of the first 
model. 

We compare this model with the first one to see that 
it provides a nice way to understand how the BPS and 
non-BPS solutions behave. We notice that 



(2) 



xx 



<t>x 



(60) 



There is only one way to make harmonic, which 

implies setting x ~~ > 0. The new first-order equations 



correspond to the equations (|23|) and 
the limit \ — ► 0. They give 



I for W {2 



i) 



This is solved by 



(x) = — tanh(Aa;) 



(61) 



(62) 



The pair of solutions formed by % = and cj> as in Eq. <p2\) 
is the pair of non-BPS solutions that we have already 
found for this model. This is another example where one 
gets non-BPS solutions as solutions to first-order differ- 
ential equations. This is perhaps the most natural way 
of understanding why this non-BPS solution has energy 
evenly distributed into its gradient and potential por- 
tions. In terms of energy densities we have 



= P2(x) = iA 2 sech 4 (Aa;) 



(63) 



and define a BPS sector, 
given by 



together with \ 
tions to get 



The explicit BPS solution is 



(68) 



0. We deal with the first-order equa- 




X <t> 



A 



rX 







(69) 



This restriction constrains the orbits in the (jp, \) plane 
in such a way that no finite orbit can connect the two 
minima, unless in the case \ — 0, in which one repro- 
duces the topological solution just obtained. This model 
has a single topological sector, and presents no family of 
topological solutions, like the family (42) and ( f43| ) ob- 
tained in the first model. 
We use to obtain 

W$=2tf>, Wg>=-2vut>, W<® = -2nx (70) 

We see that is harmonic. This means that we have 
another pair of first-order equations that come from (|22 
and (|2l). It is 



d<j> 
dx 
dx 
dx 



2fj,(j>x 



-A - fi\ + v4> 



(71) 
(72) 



We think of this pair as appearing from another function, 
as defined in the Eqs. (p6h and @. We get 



w®> = -x x -\ 



MX 



M0 2 X 



(73) 



C. The third model 

Let us now consider another model, defined by the su- 
perpotential 



This model has the first-order equations 

d( P , . ,2 2 
= -A + fi<p - fix 



dx 
dx 
dx 



(64) 

(65) 
(66) 



The model has only two minima, and for X/fi > they 
are (±\/X/fi, 0). They give 



(67) 



We notice that vanishes along the line x = where 
the former sector was explicitly constructed. 

We can use to write W$ = 2fi X , W^> = -2fix, 

and wffc = 2fi(j), showing that is also harmonic, as 
expected. As we have suggested in Sec. |^ we use the 
complex superpotential — + iW^ and the 

complex field tp = <f> + ix- We change notation and write 



Wafo) = Xtp--fi<p 3 



(74) 



as the superpotential of the third model. The new no- 
tation uses the subscript N in Wjv to identify the com- 
plex superpotential that presents the Zn symmetry - see 
Sec. VB for further details. We use W2 to get, in the vac- 
uum states (f>k,k — 1,2 



W 2 (^ 2 ) 




(75) 
(76) 
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or better 



The minima are given by (4>i, Xi), * — 1,2,3, where 



W 2 (^ fe ) = -^A^e 2 "V fc=l,2 (77) 

We use this result to obtain the tensions 

tf = \W 2 {^)-m{(Pk)\ (78) 

In the present case there is a single topological sector, 
with tension 



(79) 



Here we use tn ^ to identify tensions in models with 
com plex superpotentials having the Zn symmetry - see 
Sec. VB for more details. 



D. The fourth model 

We consider another model, defined by the superpo- 
tential 

W<*> (0, X ) = X4> + fix - \^ - Ja X 4 + \\4> 2 X 2 (80) 

This model is similar to models already investigated in 
Ref. . The equations of motion are 



= 6/xA <bx ~ 3AA(0 2 - x 2 ) - 6A 2 0( X 2 - 30 2 ) X 2 



+3A' 



2i/i2 q,.2wj,2 2 



3xW~X 2 ) 



(81) 



= 6 AA0x - 3/lA( X 2 - </> 2 ) - 6A 2 2 (0 2 - 3 X 2 )x 



+3A 2 ( X 2 -30 2 )(x 2 -0 2 )x 



(82) 



These equations are solved by field configurations that 
obey the first-order equations 



A - X^ 2 - 3 X 2 ) 
M - Ax(x 2 - 30 2 ) 
The minima of the potential obey 



dx 
dx 
dx 



0(0 2 - 3 X 2 



A 
A 



x(x 2 -W 2 ) = f - x 



(83) 
(84) 

(85) 
(86) 



Interesting cases are obtained with p = 0, and with A = 
0. They give similar models, and we consider the case 
p = 0. We use 4> 3 = A/A, and the minima now obey 



a -ix 2 ) = F 

X (x 2 ~ 3^ 2 ) = 



(87) 
(88) 



1 



Xi =o 

:<f>, X2 = ~Xa 



V3- 



(89) 
(90) 



They form a triangle in the {(f), x) plane. The ratio A/A 
specifies the vacuum states. For instance, for A = A we 
get the minima at 



(1,0) ; 



(91) 



The minima are equidistant from each other. The 
triangle they form is equilateral, and the symmetry is 
the Z% symmetry. The distance between the minima is 
d = y/3\4>\, and depend on the ratio A/A. The function 
is changed to, for p = 0, 



= A> - ^A0 4 - ^A X 4 
4 4 



2,, 2 



(92) 



The equations of motion are solved by the first-order dif- 
ferential equations 



dx 
dx 
dx 



= \[4> 3 - <j>{4> 2 - 3 X 2 )] 

= -Ax(x 2 - 30 2 ) 



We use Eq. (|9J) to get 

^ (5) (^2,X2) = W^H^xs) = -|A0 4 



(93) 
(94) 

(95) 
(96) 



We notice that W {4) (and also W^) is har monic. This 
means that the potential of this model can be obtained 
by another superpotential, given by 



(0, X ) = -H + Ax - A0 3 x + A0 X 3 
We set p = to get 

W^{4>,X) = A X - A0 3 x + A0 X 3 
The corresponding first-order equations are 



^ = Ax(x 2 -30 2 ) 
dx 



(97) 



^ = A[0 3 -0(0 2 -3x 2 )] 



(99) 
(100) 



Evidently, the minima are ((pi, yj), i = 1,2,3, the same 
minima already obtained in Eqs. ( p9|) and (|90|). However, 
behaves differently 

^ (5) (0i,Xi) = (101) 
H ?(5) (02,X2) = -H ?(5) (03,X 3 )= |V3A^ 4 (102) 
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This model leads to the model recently considered in 
Refs. We proceed like in the form model to get 

the complex superpotential 



w 3 M = \ip- \x v A 



(103) 



The field ip is complex, ip = <f> + i\. In Rcf. one 
considers the case A = A = 1. In this case the potential 
is a particular case of the potential introduced above. 
The superpotential can be seen as 
and this allows obtaining 



w 3 (?i) = - A M 4 



(104) 
(105) 

(106) 



or better 



W 3 (ip k ) = -X^e^ i ^-, k — 1,2,3 (107) 
We use this result to write 

tfk = IWsfe) - m(Vk)\ (108) 
The three tensions degenerate to the single value 



,(3) 



r\/3|A|0 4 



(109) 



IV. STABILITY 

In this Sec. [j^ we deal with issues concerning stability 
of the topological solutions in (1, 1) dimensions. We fol- 
low Ref. |f|] , recalling that the linear stability also shows 
how the bosonic fields behave in the background of the 
classical solutions. We examine linear stability for two 
coupled real scalar fields using cj>{x, t) = 4>{x) + r](x, t) 
and x( x , t) = x( x ) + £( x > Here 



r)(x, t) = 22 r) n (x) cos (w n t) 

n 



(110) 

an) 



are the fluctuations around the static solutions <f> and x- 
We use the equations of motion to get the equation for 
stability 



-1 



dx 2 



U 



(112) 



For w"^ > (w^ < 0) we have stable (unstable) solutions. 
Here 1 is the 2x2 identity matrix and 



U = 



4>4> 

Vvv 



(113) 



The bar over the potential means that after obtaining 
the derivatives, we should substitute the fields <f> and \ 
by their classical static values <f>(x) and x(x). 

We consider the first model first. In this case the ex- 
plicit calculation gives to the matrix E/W the elements 
Ujj = l^rdij, where 



TOll 



TO22 



A 



m 2 i 
A 



A 
/'■ 

- i -+2 : + >! \ 



We use the first-order equations 
duce the first-order operators 



ax 



2[i 



- + 2W (H4) 
M / 

(115) 
(116) 

|) and (|6|) to intro- 
(117) 



and 



St 



dx 



2fj, 



7,9 X 



A 
/< 

X 



(118) 



These operators allow writing Hi — S^S, where Hi is 
the Hamiltonian corresponding to the first model 



(119) 



This factorization Hi — S^S shows that Hi is positive 
semi-definite, and so can have no negative eigenvalue. We 
then conclude that the BPS solutions are stable, since 
they solve the first-o rder equat ions we need to get the 
first-order operators (117) and (118). 

In this first model, the one-component solutions can 
be solved explicitly. For the BPS pair @) and Q) we 
get the following set of discrete eigenvalues: in the first 
direction, the 6 direction 



,(1,1) _ 



= 



(1,1) 



= v3A 2 



and in the second direction, the x direction 



(i. 2 ) n 
Wq ' = 



For the pair of non-BPS solutions given by 
in Eq. ( fl5[ ) we get 



(120) 

(121) 
= and x 




2 \ 1/2 



(122) 
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and 



r,(1.2) 







,(1,2) _ 



(123) 



Since A//i > 0, we find the restriction < A//J < 1 
on the ratio A//z in order to make the non-BPS solution 
stable. 

We notice that for /i = A the eigenvalues w„ change 

I (1 2) 

to and v3A 2 , and degenerate to the values ?«„ ,n — 
0, 1, given above. This means that the scalar matter 
splits into two equal portions that present equal BPS 
states. 

Let us now consider the second model. For the pair 
of non-BPS solutions given by Eqs. ( |57| ) and (5^) we get 
the following set of eigenvalues for the bound states, in 
the first direction 



This pair of solutions represents a straight line or bit i n 
the (0, x) plane, connecting the two vacua (±1, — ^/l/3). 
We consider stability of the above solutions. Here the 
fluctuations decouple and we obtain 



U n = 9A 2 [4-6sech 2 (3Ax)] 



U 22 = 9A 2 



1 + tanh(3Aa;) sech 2 (3Aa;) 



(132) 
(133) 



The first potential U\\ is the modified Posch- Teller po- 
tential. It presents two bound states, with energies 



wo = 0, 



(134) 



The second potential U22 is not reflectionless. It supports 
no bound state, and the continuum starts at zero energy. 



.,(2,1) 



0, 



,(2,1) 



In the second direction we have 



,,(2,2) _ 



= V[(%/17-1)/2]A 2 
f' 2) = V /[(3yi7-5)/2]A 2 



(124) 

(125) 
(126) 



Let us now consider the third model. We use the pair 
obtained with (|6^) and x = to obtain for the bound 
states the set of eigenvalues, in the first direction 



it', 



(3,1) 




0. 



(3,1) 



and, in the second direction 

.,( 3 < 2 ) - / 



w. 







(127) 



(128) 



We cannot present explicit investigations concerning 
stability of solutions of the fourth model, because we do 
not know any explicit solution in this case. However, 
we consider stability of the solutions of another model, 
introduced in Ref. |24|| . This model is defined by the 
superpotential 




(129) 



The system presents the three minima (0,±y / 4/3) and 
(±1, — -\/l/3). They give three BPS sectors, having ener- 
gies |A|, 3 I A I , and 4 |A|. The highest energy is associated 
to the sector connecting the minima (±1, — yl/3). This 
is the only sector where we can find explicit solutions. 
They are given by 



4>(x) = — tanh(3Ax) 
x ( x ) = - y /T/3 



(130) 
(131) 



V. FURTHER CONSIDERATIONS 

The models investigated above present similar BPS 
and non-BPS solutions. They provide concrete ways of 
investigating how the scalar fields behave in the back- 
ground of BPS and non-BP S solu tions. We ex plore some 
possibilities in the next Sec. V A. In Sec. |VB| we investi- 
gate intersection of defects. 



A. Entrapment of the other field 

In the first model, we can further understand the im- 
portance of the topological solutions by investigating the 
potential V\ ((f), x) ■ It is given by 



V 1 (^x) = \>?(4> 2 



l) 2 + AM0 2 - 1)X 2 



(135) 



This model presents the one-component BPS solution 
© and @. We have 4> 2 = 1 outside the BPS de- 
fect, and 4> = at x — 0, at the center of the defect. We 
then use V\ (</>, x) to get 



V-i 



2 A 

x — 



Lx) 



1 2 4 
2« 



(136) 
(137) 



For A//i > the potential (136) shows the presence of 
spontaneous symmetry breaking in the independent x 
direction in the (0, x) plane. In this case the x 
presents squared masses given by 



m 2 (in) = AX/i 
m 2 (out) = A^i 2 



(138) 
(139) 



where in and out identify the region inside and outside 
the defect. We see that the ratio of masses is 







in; 



.(in) 



'{out) 



X 
/' 



(140) 



Here we have 



This ratio identifies the region 0<A//i<lin the space 
of parameters where the field 4> entraps the other field. 
For a given A, we see that the parameter /i controls the 
ratio of squared masses. Thus, the efficiency of the en- 
trapment depends on u. 

This model can be used as a bag model |^9fl , similar 
to the models in Refs. pO| , pl| . Here we can further en- 
large the model by changing dimensions from (1,1) to 
(3, 1) and making the x field complex, going from the 
Z2 x Z2 symmetry to the Z2 x U(l). However, since in- 
side the defect formed by (j)(x), the x field still engenders 
spontaneous symmetry breaking for Xu > 0, we can only 
circumvent the presence of Goldstone bosons when one 
gauge the related t/(l) symmetry, making it local. The 
model that emerges presents the symmetry Z2 x Ui(l). 
Since we started in (3, 1) dimensions, inside the wall the 
theory is effectively (2, 1) dimensional. This is an al- 
ternate mechanism for dimensional reduction, and inside 
the wall the effective theory may very naturally develop 
the Callan-Harvey effect [Q, if one appropriately cou- 
ples fermions to the second field, x- This scenario seems 
to appear interesting, and we shall further explore this 
possibility in a separate report. 

The first model is simpler for Xu < 0. In this case, 
spontaneous symmetry breaking no longer appear inside 
the wall formed by the (j> field. However, the wall is still 
of the BPS type. Th e po tential V\(4>, x) changes in a way 
such that the ratio ( |l40| ) becomes 



m: 



\ (in) 



'■(out) 



X_ 

2jl 



(141) 



This result shows that the topological defect generated 
by (j) can still entrap x mesons for > A//i > —2, acting 
like a bag for the elementary excitations of the field x- If 
we change the symmetry Z2 x Z2 to Z2 x U(l) by mak- 
ing the x field complex, in this case we do not need to 
make the U(l) symmetry local, and we can still explore 
the presence of a global symmetry within the wall. Work 
in this direction was already presented in Ref. |41|, in 
a model inspired in the work of Ref. Q, with a poten- 
tial that differs from the one here proposed. We believe 
that our model may give further insight into the issues 
that arise in this scenario. Further applications can be 
done, as for instance in Ref. p2| ], where one explores the 
presence of diffuse domain walls at intersecting flat di- 
rections within the cosmological scenario. Also, we can 
follow Ref. jl3| to investigate issues concerning hybrid in- 
flation, envolving features that requires at least two real 
scalar fields, the inflaton field and another one, which is 
in general used to break the symmetry. 

The second model is defined with the superpotential 
([52|). It gives the potential 

VK0, X ) = iA 2 (0 2 -l) 2 + 2AW (142) 



1 



v 2 (4> = o, x ) = 2 a 2 



(143) 

^ 2 (0 2 = 1,X) = 2AV (144) 
Thus the squared masses associated to the x field are 

m 2 (in) = (145) 



m 2 (out) = 4A 2 



(146) 



In this case the non-BPS defect always entraps the field 
X, acting like a bag for the elementary x excitations. 
Here, however, the x meson is massless inside the bag. 
This shows that the entrapment is more efficient in this 
model, although the basic solution, the bag, is of the 
non-BPS type. 

We can investigate the third model similarly. It is de- 
fined by the superpotential ( |64| ) and we have 



1 



- V - 



1 2 4 1 n 2 i2 2 
-U X +2/1 (j> X 



We can write, for X/u > 0, 
1 



2 V M 



2, A 



V 3 (±VX/u,x)^2X^ X Z + ^^X 



(147) 

(148) 
(149) 



No spontaneous symmetry breaking appears inside the 
</>-defect anymore, and so the model cannot support de- 
fect inside defect. However, the squared masses of the 
elementary x mesons are 



m 2 (in) 



2A/i 



fh 2 (out) = 4A/i 



We get the ratio 



~ 2 
X(out) 



(150) 
(151) 



(152) 



This ratio does not depend on the parameters A and /z, 
but the topological defect may also entrap x mesons. 
This is another model, in which the topological state may 
simulate a bag to entrap the elementary excitations of the 
other field. In this case, however, the efficiency of the en- 
trapment can no longer be controlled by the parameters 
that defines the model. 

The second and third models may also be immersed 
in the (3, 1) dimensional space-time. This immersion 
provide alternative ways to investigate the behavior of 
the current generated by making the second field com- 
plex, as commented above. They give similar scenarios 
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for charged walls, and may present distinct physical con- 
tents. 

Other applications include the diverse issues recently 
investigated in Refs. [^2|- ^5| , p7p^ 1 . We can, for instance, 
enlarge the models in order to break the symmetry 
and/or the supersymmetry. This possibility may give rise 
to Fermi balls [^7), Fermi disks (TJJ], and other charged 
objects |44|l , which depend on the particular space-time 
dimension one starts with, and on the parameters that 
define the terms added to allow the symmetry breaking. 
We can also follow Refs. |l2|,[l3|,[l5| to extend the scenario 
of nested defects to the case of nested network of defects. 
This picture follows as a direct generalization of the re- 
sult of Ref . , which shows explicitly that defects may 
nest defects if the model engenders the Z2 x Z2 symmetry 
in systems of two real scalar fields. As we are going to 
see below, to have intersection of defects we must neces- 
sarily enlarge the Zi symmetry to at least the Z% one. 
For this reason, we can think of a model that presents 
the Z2 x Z% symmetry, composed of three real scalar 
fields, one presenting the Z 2 symmetry, and the other 
two controlling the Z 3 symmetry. In (3, 1) space-time di- 
mensions the first field engenders the Zi symmetry, and 
we can certainly get to the case where it may give rise 
to a domain wall, which may nest the two other fields 
in its interior. Within this scenario, if the Z3 symmetry 
related to the second and third fields is still effective we 
may nest a planar network of defects inside domain walls. 
We take advantage of the model investigated in Ref. ]24| ] 
to present a model that seems to work correctly. This 
model is defined by 



V M x) = A 2 (0 2 + xT - x<KF - 3x 2 ) 

-\\\4> 2 +X 2 ) (157) 

We see that the Z3 symmetry of the pair of fields (4>, x) is 
preserved both inside and outside the domain wall gener- 
ated by the first field, tp. In the region inside the domain 
wall, Vi n {4>, x) has minima at the points 



(158) 



In this region inside the domain wall the squared masses 
corresponding to the two mesons are 



27 

m l( in ) = -J* 2 



(159) 
(160) 



In the region outside the domain wall, V ou t (4>, x) has the 
minima [E4fl 



.0 



14* 



(161) 



In this case the squared masses are m^out) — nr^(out) = 
(27/2) A 2 . These results show that the field ip engenders 
an efficient mechanism for the entrapment of the pair of 
fields (4>, x) m this model. We shall further explore this 
and other related issues in another work. 



F(^,0,X)^A 2 U 2 ~\ 



aV (0 2 + r 

+A 2 (0 2 + x 2 ) 2 -A 2 0(0 2 -3x 2 ) (153) 



It behaves standardly in (3, 1) space-time dimensions and 
is such that V(tp) = V(ip,0,0) obeys 



5 V I * 



9 



(154) 



We consider that the field ip = ip(x) depends only on x. 
Thus, it supports BPS defects with the explicit form 



tanh 



(155) 



The corresponding energy density is 3V3A 2 . We also 
notice that the potentials V* n (</>, x) = V(Q,<j),x) and 
V ou t{<t>,x) = V(±3/2, <j>,x) are such that 
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,X) = | A 2 + A 2 (0 2 + X 2 ) 2 -AW-3x 2 ) 

(156) 



and 



B. Intersection of defects 

The process of intersection of defects was already ex- 
plored in Ref. p0| , in the case of supersymmetric field 
theories. There, the basic idea was to show that the 
intersection of two extended objects is energetically ad- 
missible when it can be seen as a reaction that occurs 
cxothcrmically. This may be translated into the expres- 
sion tik < Uj + tjk , where Uk , Uj , and tjk are the ener- 
gies or tensions associated to the fusion of the two de- 
fects, and to the individuals defects, respectively. The 
crucial point here is that in supersymmetric theories the 
topological charge T of the defect appears as a central 
charge, and for BPS states one gets t = |T|. If one uses 
hj = \%j\,tjk = \Tjk\, and t ik = |7y +T jk \ the triangle 
inequality Uk < Uj + tjk follows naturally. This triangle 
inequality is strictly valid in supersymmetric systems de- 
scribed by complex superpotentials, whose values at the 
vacuum states define points in the complex plane, not 
aligned into a straight-line segment. 

The fourth model is the supersymmetric model that 
opens the possibility for the presence of stable intersec- 
tion of extended objects, giving rise to a triple junction 
of BPS states, as recently shown in Ref. ]2lJ - see also 
Refs. [E2lE3l. This model was considered in Sec HID, and 
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the basic model is given by ip — (1/4) ip A . Here, however, 
we consider the case 



1 



N + l 



(162) 



with A and A real parameters. This model presents the 
Zn symmetry. The vacuum states are singular points of 
the superpotential. They obey ip N — A/A, and the N 
solutions are, using A = a N A with a real and positive, 



a e 



2-Ki[(k-l)/N] 



1,2, 



,N. 



(163) 



The topological sectors connect pairs of adjacent vacua, 
and in general the number of BPS sectors is N(N — 1) /2. 
Their energies or tensions are given in terms of the values 
of the superpotential in the singular points, that is 



t 



(JV) 



|VM^)-VM^-)I 



(164) 



The tensions can be classified according to the integer n, 
n = 1, 2, < [N/2], where [N/2] stands for the greatest 
integer not greater than N/2 itself. These tensions have 
the explicit values 



tW = — lAI 

iV + i |A| 



1 — cos ( 2n 



N 



(165) 



The integer n identifies topological sectors described via 
the connection between a given vacuum, and its n — th 
neighbor: n — 1 for connections between first neighbors, 
n = 2 for second neighbors, and so forth. We can check 
that these tensions are ordered in the form < 4;+!) 
for k = 1, < [N/2], and N > 4. This order also shows 
that the direct connection between two vacuum states is 
always less energetic than any other possible connection. 



For instance, t { 2 N) < 2t\" '' , TV > 4, and ^ Jv; < t 



(N) 



AN) 



t[ N) < 3t\">, N > 6, and so forth. 

More recently, in Rcf. |24|] we considered another 
model. The model is defined by the potential 



AN) 



2\ 2 



27, 



■tf-\*4>W-Ztf) + —X' (166) 



It presents the Z% symmetry, but it is described by a 
potential that cannot be written in terms of some super- 
potential. In this case the key issue relies on finding a 
way out of supersymmetry to show that the fusion of de- 
fects still occurs exothermically. In the supersymmctric 
case one can use interesting properties to make the rea- 
soning naturally simple. In the non-supersymmetric case, 
however, we have found non-BPS solutions that presents 
the interesting property of having gradient and poten- 
tial portions of the energy evenly distributed to compose 
the total energy of the defect. As we have shown, such 
a property is also shared by all the BPS solutions, and 
so we can use it to present a new reasoning, that works 



within and beyond the context of supersymmetry. The 
reasoning follows from the topological current presented 
in Eq. (pOL). Here we have 



p(x) 



<t>'{x) 



(167) 



Thus /?' p = 2g(x), and for static solutions that presents 
equal gradient and potential portions we get p 1 p — g(x) + 
p(x) = e{x), or better 



dx plj pij 



(168) 



We consider fusion of two defects in the model investi- 
gated in Ref. . The topological defects there obtained 
obey g(x) = p(x). Thus, we can use 



(pij + Pjkfipij + Pjk) = P%Pij 



PjkPjk 



+PijPjk + PjkPij 



(169) 



Also, we take advantage of the symmetry of the model to 
see that if the pair (0, x) represents a solution in a given 
sector, all the other solutions can be obtained by rotating 
the solution according to the Z% symmetry. This means 
that we can write 



cos(oiy) sin(aiy) 
-sin(aij) cos(ary) 



(170) 



Here ctij is the angle between the given solution (4>,x) 
and the solution (<fiij,Xij)- I n ^ ne ^3 model under con- 
sider ati on a g can only be 27r/3 and 47r/3. We now use 



(168), (169), and (|170|) to obtain t lk < t y +t jk . This 



result shows that the three-junction in the model under 
consideration occurs exothermically. 



VI. NETWORKS OF DEFECTS 



The recent Refs. |21|-|23fl have introduced investigations 
concerning defect junctions and networks in supersym- 
metric systems, and now we turn our attention to this 
possibility. We consider junctions in the plane, which 
may give rise to planar networks of defects. To imple- 
ment this possibility we work in the (2, 1) space-time, in 
the plane (x, y). We identify the axes (x, y) with the axes 
(<j), x) of the space of configurations. We illustrate this 
situation b y con sidering the model defined by the poten- 
tial in Eq. 



case arc 124 



166 ). The topological solutions in the planar 



6 (2,3) 



X%,3) 



3 

_ 4 
3 

± — \/3 tanh 

4 




(171) 
(172) 



and 



12 



X(l, 2 ) 



1 T 3 tanh j i W y A (y + V3x) 



;v5 



1 ± tanh ( i J y A (y + v^x) 



and 



1 ± 3tanh | A (y - \/3x) 



X 



(1,3) 



V3 



1 =F tanh 




V3a 



(173) 



(174) 



(175) 



(176) 



They may be used to represent the three-junction in the 
thin wall approximation. 

Let us consider the first model, in the case A//i > 0. 
The minima are at v\ = (1,0), V2 — (0, \J X/fi), = 
(—1,0), and W4 = (0,—y/X/fj,). The energies or tensions 



of the BPS and non-BPS sectors are t BPS 



2t 



BPS — 

BPS — (4/ 3) | A | . For a given A, the value of 
fi allows obtaining three distinct situations. The first 
concerns the case /i > A. In this case the minima give 
rise to an array of defects of the form shown in FIG. 1, 
in the thin wall approximation. 



network in the two following ways. The s- network, which 
appears for X/fi < 1, and the t-network, for X/fi > 1. 
They are depicted in FIG. 4 and in FIG. 5, in the case of 
regular hexagonal networks in the thin wall approxima- 
tion. In the first case we impose that the triangle with 
vertices Vi,V2,V4 is equilateral. This gives A//i = 1/3. 
In the second case we impose that the triangle with ver- 
tices vi,V2,va is equilateral. It gives A//i = 3, the result 
we should expect from duality. We notice that the first 
(second) imposition subdivides the distance between the 
horizontal (vertical) vacua into three equal pieces, which 
is the condition for the formation of a network of regular 
hexagons. 




FIG. 2. A nother array of defects, suggested by the first 
model of Sec. 



[II A in the case of X/^l = 1. 





FIG. 1. A possi ble a rray of defects, as suggested by the 
first model of Sec. [II A in the case of X/fi < 1. The vacuum 
states are represented by vi = 1,V2 = 2,vg = 3, and V4 = 4, 
and the thinner and thicker lines stand for BPS and non-BPS 
defects, respectively. 

The second case concerns /i = A, and the third one 
fi < X. In the thin wall approximation they are shown 
in FIG. 2, and in FIG. 3, respectively. The second case 
seems to give an intersection of four half-defects, but this 
is not true since for /i = A the two fields decouple. Thus, 
the apparent intersection of four half-defects represents 
in fact the uninteresting case of two non-intersecting de- 
fects. 

The first and third cases seem to be dual to each other, 
since they are linked by the interchange A *-* /x. This 
duality is very much like the s <-> t duality that the dia- 
grams depicted in FIG. 1 and in FIG. 3 remember. This 
scenario excludes the self-dual case, where /i = A. 

The basic diagrams of FIG. 1 and FIG. 3 can be used to 
generate planar networks of defects. We can generate a 





FIG. 3. Another po ssible array of defects, as suggested by 
the first model of Sec. 



[II A 



in the case of X/n > 1 . 



For \x = 3A, we have two sides of the hexagon formed 
by non-BPS states, with tension \/{l/3)t BPS . The 
other four sides are formed by BPS states, with tension 
{l/2)t 1 BPS . For /i = (1/3)A, two sides of the hexagon are 
formed by BPS states, with tension t BPS , and the other 
four sides are also formed by BPS states, with tension 
(1/2) t BPS . 

We think of a planar network of defects. We notice 
that a t-network means X/fi > 1, and this allows hexag- 
onal networks with two distinct BPS states, having two 
different tensions, t BPS and i BPS = (l/2)t BPS , show- 
ing that each three-junction of this t-network is at the 
threshold of stability. For A/jU < 1, in the case of s- 
networks, however, the situation is different since now 
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we must have two non-BPS states to build the hexago- 
nal cell. Since the ratio A//x < 1 controls the energy of 
the non-BPS state, every three-junction may obey the 
inequality U 3 < t jk + t ki . 




FIG. 4. The s-network, depicted as a regular hexagonal 
network of defects in the case of X/fi — 1/3. The thicker lines 
represent non-BPS states. 




FIG. 5. The t-network, depicted as a regular hexagonal 
network of defects in the case of A//i = 3. 



We can consider the case where t 



nBPS 



= t 



BPS- 



This 



gives /i = 4 A and the tensions degenerate to the single 
value (2/3)|A|. The hexagonal cell is no longer regular. 
The angle between the two BPS states changes from a — 
120 to a = 2 9, where 9 is slightly lesser than 63.5 degrees. 
We notice that we cannot make the tensions and angles 
degenerate simultaneously, and this indicates instability 
of the suggested networks. 

Effects that contribute to stabilize the network may 
appear when we go beyond the classical level and intro- 
duce the quantum corrections, since the quantum con- 
tributions may modify the classical scenario. For in- 
stance, in the case A//i > 1, in a t-network all the defects 
arc BPS defects, and their energies or tensions are pro- 
tected against receiving quantum corrections. On the 



other hand, for A//i < 1, in the case of a s-network there 
are non-BPS states whose tensions are not protected and 
may receive quantum corrections, changing the classical 
value t^gpg, contributing to stabilize the network. 

The first model does not allow the appearence of a net- 
work of planar squares, although squares can also fill the 
plane. In our model this would require the symmetry, 
which is not effective in that model. 

The fourth model leads to the model considered in 
Refs. pljjffi] . In this case the three-junction is formed by 
BPS defects, and we can form an hexagonal array. Here, 
the three-junction is stable, because each leg carries the 
same tension. But it does not generate a BPS network, 
since any two adjacent three-junctions must have differ- 
ent winding to accomodate the three vacuum states in 
the network, and this would imply that each side of the 
hexagonal cell must be seen as BPS and anti-BPS, simul- 
taneously, frustrating and destabilizing the network - see 
Refs. |lj-||. 




FIG. 6. The regular hexagonal network of non-BPS d efect s 
given by the Z$ model defined by the potential of Eq. (166). 



We consider another model, investigated in Ref. pi[ . 
This model is defined by no superpotential, and the argu- 
ment above is no longer valid. The Z3 symmetry seems 
to allow the presence of a network of defects, also in the 
form of regular hexagonal array. In the thin wall approx- 
imation, the hexagonal network that appears in this case 
is interesting because the solutions are explicitly known, 
form straight-line segments joining the vacuum states, 
and have the same energy or tension, (9/4)^/27/8 |A|. 
The Z3 symmetry that appears in the vacuum sectors 
is also effective in the sectors of the topological defects. 
The inequality f jj. < tfj + tj k is strictly valid in this case, 
ensuring stability of the hexagonal network, as depicted 
in FIG. 6 in the thin wall approximation. 



VII. COMMENTS AND CONCLUSIONS 

In the present work we have examined several mod- 
els of coupled scalar fields. We have found some explicit 
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solutions, of the BPS and non-BPS type. We have also 
investigated linear stability, checking the behavior of the 
bosonic matter on the background of the BPS and non- 
BPS defects. The BPS states are solutions that obey 
first-order differential equations, and the non-BPS states 
are not. The BPS states engender the property of having 
energy evenly distributed in the gradient and potential 
portions. However, in the systems of two coupled fields 
that we have introduced, the non-BPS states also engen- 
der the same feature, as we have explicitly shown in this 
work. The point here is that the non-BPS states of the 
systems of coupled fields are BPS states of simpler sys- 
tems, which are projections of systems of coupled fields 
onto a single, specific field direction. 



The first model is the model investigated in Sec. Ill A 



It presents several properties. In particular, it has a rich 
BPS sector, the sector that connect the minima (1,0) 
and (-1,0) by one-component and two-component BPS 
states. The two-component state describe a family of so- 
lutions, parametrized by the ratio between the coupling 
constants that specify the system. These one-component 
and two-component states allow diverse applications, like 
the applications as bag models in Refs. p7| ], as defects in- 
side defects in Refs. [ pd|Jl^ , p^ , ^5| , and in applications to 
condensed matter, to polimeric chains having two rel- 
evant degrees of freedom to represent the chain - see 
Refs. [|3|. 

In applications in the form of bag models, in addi- 
tion to the models already introduced in Refs. |3(],|3l| we 
have other possibilities, given by the first, second and 
third models introduced in Sec. III. These models dif- 
fer from the model considered in Ref. jy], and we be- 
lieve that this new direction may shed further light on 
the issues explored in that work. In particular, it seems 
that our suggestions are easier to examine, and may offer 
an analytical understanding of the issues put forward in 
Ref. J4l| under numerical investigations. We recall that 
the model investigated in jll| is defined by a potential 
that cannot be reduced to none of the models considered 
in the present work. 

In application to defect inside defect, much work was 
already done in Refs. [[fl 15 1. The natural extension is 
now related to issues where one allows the presence of lo- 
cal symmetries, as in the work on superconducting strings 
||. In the first model of Sec III the Z2 x Z2 symmetry 
can be enlarged by enlarging each one of the two inde- 
pende nt sy mmetries. As we have already commented on 
in Sec. 



Ill 



we can for instance make the \ field complex, 
in this case changing the Z2 x Z2 symmetry to Z2 x Ui(l), 
if one also includes the Abelian gauge field. This opens 
new possibilities, providing for instance a natural sce- 
nario to the Callan-Harvey effect Another possibil- 
ity includes models that presents the Z 2 x Z 3 symmetry, 
which may provide scenarios for domain walls that nest 
planar networks of topological defects. These issues are 
presently under consideration, and we hope to report on 
them in the near future. 

We have also commented on junctions of defects, and 



the subsequent formation of networks of defects. In this 
case we need to enlarge the spatial dimension from d = 1 
to at least the d = 2 case. In the planar case, we have 
exemplifyed how to get planar solutions from the unidi- 
mensional kinks of the model. Here we have considered 



three models, the first and the fourth model of Sec. Ill, 
and another one, firstly considered in Ref. [ p4j . These 
systems are all different: the first model admits a super- 
symmetric extension of the type N = 1, and seems to 
admit no stable network of defects, at least classically; 
the fourth model presents junctions that break 1/4 su- 
persymmetry, but it also seems to admit no stable net- 
work of BPS defects. The other model does not present 
a supersymmetric extension, but the defect solutions get 
even contributions from the gradient and potential por- 
tions of the energy, and so they behave very much like 
BPS states. This is crucial for showing that the fusion of 
defects is a reaction that occurs exothcrmically, and so 
the three-junction generates a stable network of defects. 
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